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Quark Loop Contributions to Neutron, Deuteron, and Mercury EDMs from 

Supersymmetry without R parity 



We present a detailed analysis together with numerical calculations on one-loop contributions to 
the neutron, deuteron and mercury electric dipole moment from supersymmetry without R parity, 
focusing on the quark-scalar loop contributions. Being proportional to top Yukawa and top mass, 
such contributions are often large, and since these are proportional to hitherto unconstrained com- 
binations of bilinear and trilinear RPV parameters, they are all the more interesting. Complete 
formulas are given for the various contributions through the quark dipole operators including the 
contribution from color dipole operator. The contribution from color dipole operator is found to be 
similar order in magnitude when compared to the electric dipole operator and should be included 
in any consistent analysis. Analytical expressions illustrating the explicit role of the R-parity vio- 
lating parameters are given following perturbative diagonalization of mass-squared matrices for the 
scalars. Dominant contributions come from the combinations B* X'iji for which we obtain robust 
bounds. It turns out that neutron and deuteron EDMs receive much stronger contributions than 
mercury EDM and any null result at the future deuteron EDM experiment or Los Alamos neutron 
EDM experiment can lead to extra-ordinary constraints on RPV parameter space. Even if R-parity 
violating couplings are real, CKM phase does induce RPV contribution and for some cases such a 
contribution is as strong as contribution from phases in the R-parity violating couplings. Hence, 
we have bounds directly on |_B*Aiji| even if the RPV parameters are all real. Interestingly, even 
if slepton mass and/or no is as high as 1 TeV, it still leads to neutron EDM that is an order of 
magnitude larger than the sensitivity at Los Alamos experiment. Since the results are not much 
sensitive to tan /3, our constraints will survive even if other observables tighten the constraints on 
tan /3. 

PACS numbers: .. 



Chan-Chi Chiou and Otto C. W. Kon^!! 



Department of Physics, National Central University, Chung-Li, TAIWAN 32054- 



Rishikesh D. Vaii 




Department of Theoretical Physics 
Tata Institute of Fundamental Research, 
Mumhai 400005, India 



* Electronic address: 'otto@phy.ncu.edu.tw| 
t Electronic address: |rishidilip@gmail.com| 



2 



I. INTRODUCTION 



The problems of neutrino mass, baryogenesis, dark matter, dark energy and gauge hierarchy, provide unambiguous 
hints toward physics beyond standard model (SM). Whereas a direct discovery of new physics particles at colliders 
is indispensable, a search for alternative observables could not only provide a means to discovery of new physics 
but also prove complimentary by hinting at favorable regions in parameter space. Discrete symmetries and their 
violations have been crucial to establishing and validating the SM. Forty years after the discovery of CP-violation 
its experimentally observed effects in the K and B-meson systems [l|, Q are generally compatible with the standard 
model (SM) predictions with the Kobayashi-Maskawa (KM) phase as its sole source. The search for more CP violating 
observables is keenly pursued at ongoing and upcoming B physics experiments, and it is largely confined to flavor 
changing sector. Within the flavor diagonal sector, P and T violating electric dipole moments (EDMs) of fermions 
3 , heavy atoms and molecules are interesting CP violating observables that provide essentially background free and 
sensitive probes of physics beyond SM [1,01 • Though the search for non- vanishing EDM has so far yielded null results, 
the present experimental scenario with regard to EDM measurements is very encouraging, with most of them well 
within the range of interesting predictions from physics beyond SM. For the convenience of reader, below we briefly 
describe the current status of EDM experiments. This will also serve to motivate our case for specific new physics 
contributions that we discuss here. 

Since the early work of Purcell and Ramsey [6|], EDM experiments have dramatically improved in precision. The 
standard method for measuring a permanent EDM of a particle is by placing it in an external electric field E and 
looking for a shift in energy that is linear in E. This explains the focus on the electrically neutral candidates for 
EDM measurements. Naively, the Schiff screening theorem prevents any measurement of atomic EDM. Indeed, the 
assumptions of Schiff theorem, namely the point-sized nucleus and non-relativistic limit, are significantly violated in 
heavy atoms and hence facilitate EDM measurement. Ironically, the shielding which is complete in non-relativistic 
limit, actually produces an enhancement in the realistic relativistic limit JJ. Sanders [S] pointed out that due to 
relativistic magnetic effects, the atomic EDM induced in heavy atoms can be strongly enhanced compared to electron 
EDM inducing it, leading to the best limit on electron EDM. The enhancement factor can be over two orders of 
magnitude. For the case of paramagnetic atoms, the best bound so far is for the thallium ( ^°^T1), |c?ti| < 9.0 x 10"^^ 
(90%C.L.) leading to the tightest limit on electron EDM |c?e| < 1.6 x 10^^'' ^9]. Note that the numbers are in the 
standard e cm unit, which is assumed throughout the paper. For the case of diamagnetic atoms, mercury ( ^^^Hg) 
EDM is best constrained by the Washington group, the bound being Jdngl < 2 x 10^^^ (95 % C.L.) [l3|. An upgraded 
experiment is expected to improve the accuracy by a factor of four [ll| . The violation of Schiff theorem comes about 
due to the finite size effect of the nucleus. For the nucleon EDM, the best bound so far is on the neutron EDM 
from the Grenoble experiment, |(i„| < 6.3 x 10"^^ (90% C.L.) [l^. The most recent result from the same experiment 



\dn\ < 3 X 10~2^ [13]. It is expected to reach a goal of |d„| < 1.5 x lO"^'^ [ij. The Los Alamos neutron EDM 
experiment will provide two orders of magnitude improvement, probing dn down to order 10~^* [T5| . The standard 
method of measuring the energy shift linear in electric field fails for the EDM of charged particle due to acceleration 
of charged particle. However, in recent years a new dedicated method of searching for EDM of charged particles in 
storage rings has been developed [H, [O] . A muon EDM experiment is proposed that is expected to reach a sensitivity 
of 10"^'* which is an improvement of a factor of 10^ to 10^ over the last CERN muon g-2 experiment [l^. A deuteron 
EDM experiment, again using storage ring is proposed, that would reach a sensitivity of 10~^^ that is 10 to 100 times 
better than current EDM limits in terms of sensitivity to quark EDM and QCD 9 parameter [T^. 

With so many exciting ongoing and upcoming experiments, EDM searches certainly provide complementary alter- 
native to probe physics beyond the SM. In this work we will focus on the neutron, deuteron and mercury EDM. As 
an example of physics beyond SM, we shall focus on supersymmetry (SUSY) without R-parity or the generic super- 
symmetric standard model '^I'l. When the large number of baryon or lepton number violating terms are removed 
by imposing an ad hoc discrete symmetry called R-parity, one obtains the MSSM Lagrangian. GSSM is a complete 
theory of SUSY without R-parity, where all kinds of R-parity violating (RPV) terms are admitted without bias. It is 
generally better motivated than ad hoc versions of RPV theories. The MSSM itself, without extension such as adding 
SM singlet superfields and admitting violation of lepton number, cannot accommodate neutrino mass mixings and 
hence oscillations jl^. Given SUSY, the GSSM is actually conceptually the simplest framework to accommodate the 
latter. The large number of a priori arbitrary RPV terms do make phenomenology complicated. However, the origin 
of the (pattern of) values for the couplings maybe considered to be on the same footing as that of the SM Yukawa 
couplings. For example, it has been shown in [23| that one can understand the origin, pattern and magnitude of all 
the RPV terms as a result of a spontaneously broken anomalous Abelian family symmetry model. 

As we will see in the next section, EDM of mercury, deuteron and neutron are expressible in terms of quark EDM 
and color EDM (CEDM). Neutron EDM is an old favorite and within SM, the CKM phase contribution starts at 
three loop level i24i1 . Generic ( R-parit y conserving) SUSY contributions starts at one loop level and hence can be 
large [1^, [2^, [27l. l28l. [29l . [30l . [Sll . [32 . [33l [H, [sH] ■ Most of the phenomenology studies of the case admitting R-parity 
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violation have largely been confined to the trilinear superpotential parameters. Within the latter framework, it has 
been shown that contributions to fermion EDMs start at two loop leve l 13611 . Striking one-loop contributions are, 
however, identified and discussed based on the GSSM framework l38j. Under that generic setting, all RPV 
couplings are considered without bias. Note that studies of RPV physics under such a generic setting is uncommon. 
Not only that a priori approximations were usually taken on the form of R-parity violation especially on the set of 
bilinear parameters, such approximations were often times not clearly stated, if appreciated well enough. They may 
be confused with the flavor basis choice issues (see Ref.(2l| for detailed discussions on the aspects). EDM study from 
Godbole el.al in Ref.jl^, however, did state explicitly the bilinear couplings were neglected in their study. The one 
loop contributions are exactly resulted from combinations of a bilinear and a trilinear RPV couplings. Similarly, flavor 
off-diagonal dipole moment contributing to the case of the 6 ^ 57 decay [39,] and that of the — > 67 [i^ decay at 
one-loop level have also been presented. That is the approach taken here. In Ref. [131, the one loop EDM contributions 
from the gluino loop, chargino-like loop, neutralino-like loop are studied numerically in some details. It shows that the 
RPV parameter combination /i*A'j^ dominates, with small sensitivity to the value of tan/3. The experimental bound 
on neutron EDM is used to constrain the model parameter space, especially the RPV part. However, the alternative 
one-loop contribution containing a quark and a scalar in the loop is also important. For the case of down quark 
EDM, it contains a top-quark loop and hence proportional to top mass and top Yukawa, thus can giving rise to very 
large contribution. We give complete 1-loop formulas for these contributions to the EDMs of the up- and down-sector 
quarks, with full incorporation of family mixings (in Ref . [37| . family mixing was neglected for simplicity). We present 
numerical analysis of quark-scalar loop contributions from all possible combinations of RPV parameters. Besides the 
familiar /i,*A'j^, there are a list of combinations of the type B*)Ljf. which are particularly interesting and will be the 
focus of this paper. In comparison to the earlier works [stI [38j. our investigation provides more extensive results, 
both analytically as well as numerically with elaboration on some physics issues as well as includes contributions to 
deuteron and mercury EDM. We obtain bounds on combinations of RPV couplings which are otherwise unavailable. 
In fact, as we will see, any null results of future EDM experiments lead to stringent constraints on RPV parameters 
that have not been constrained so far. 

Note that we decouple any explicit discussion of leptonic EDM, like for the electron and muon as well as c?ti 
in which the electron EDM has a dominant role, from the study here. Unlike the case of the R-parity conserving 
contributions, the RPV contributions to EDMs of the quark and lepton sectors have quite independent origin. For the 
lepton sector, it involves the A-couplings, rather than the A'-couplings. Hence, a whole set of different combinations 
of RPV parameters are to be constrained by the leptonic EDMs numbers — not to be addressed in this paper. 

The structure of the paper is as follows. In the next section we describe our notation and framework, the so called 
single- VEV parametrization (SVP). We will also describe the formula for neutron, deuteron and mercury EDM in 
terms of quark EDM and chromoelectric dipole moment (CEDM). In section HI, we will describe the quark loop 
contribution to quark EDM and CEDM coming from a combination of a bilinear and trilinear RPV couplings. In 
section IV we will discuss the results and in section V we will conclude. 



II. FORMULATION AND NOTATION 



We summarize the model here while setting the notation. Details of the formulation adopted is elaborated in 



Ref. [2l|. The most general renormalizable superpotential for the supersymmetric SM (without R-parity) can be 



written as 



+ ^al3kLaL''0Ek^ + - X-jkU^ Df Dk , (1) 

where (a, b) are SU{2) indices, (i, j, k) are the usual family (flavor) indices, and (a, [3) are extended flavor index going 
from to 3. In the limit where Xijk, A'^tj, A"^^, and /i^ all vanish, one recovers the expression for the R-parity preserving 

case, with Lq identified as Hd- Without R-parity imposed, the latter is not a priori distinguishable from the Li's. 
Note that A is antisymmetric in the first two indices, as required by the SU{2) product rules, as shown explicitly here 
with — —£21 — 1- Similarly, A" is antisymmetric in the last two indices, from SU{3)c- 

The large number of new parameters involved, however, makes the theory difficult to analyze. An optimal 
parametrization, called the single- VEV parametrization (SVP) has been advocated[23| to make the the task manage- 
able. Here, the choice of an optimal parametrization mainly concerns the 4 La fiavors. Under the SVP, flavor bases 
are chosen such that : 1) among the L^s, only Lq, bears a VEV, i.e. (Li) = 0; 2) h^jk{= ^ojk) = diag{mi, m2, ma}; 3) 

hjk{= Xojk = -^jofe) = ^d[a.g{md,ms,mb}; 4) ^^j, = ^^km diag{m„, mc, mt}, where Vg = {Lq) and w„ = y/2{Hu). 
The big advantage of here is that the (tree-level) mass matrices for all the fermions do not involve any of the trilinear 
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RPV couplings, though the approach makes no assumption on any RPV couphng inchiding even those from soft SUSY 
breaking; and all the parameters used are uniquely defined, with the exception of some possibly removable phases. 
The soft SUSY breaking part of the Lagrangian in GSSM can be written as follows : 



Q'^ml Q + U^vT^ U + D'<ml D + DmlL + E^ml E - 
+^BB + -^WW + -^99 + h.c. , 



™2 



\Hn 



h.c. 



(2) 



where we have separated the R-parity conserving ones from the RPV ones {Hd = Lq) for the A-terms. Note that 
L''rh?L, unlike the other soft mass terms, is given by a 4 x 4 matrix. Explicitly, is fh^^ of the MSSM case while 
m^^^ 's give RPV mass mixings. 

Details of the tree- level mass matrices for all fermions and scalars are summarized in Ref. |21| . For the analytical 
appreciation of many of the results, approximate expressions of all the RPV mass mixings are very useful. The 
expressions are available from perturbative diagonalization of the mass matrices [2ll|. 

III. THE QUARK LOOP CONTRIBUTION TO NEUTRON, DEUTERON AND MERCURY EDM 

As we will see toward the end of this section, neutron, deuteron, and mercury EDMs are ultimately expressible 
in terms of quark EDMs and CEDMs. Quark EDMs and CEDMs are typically defined from the following effective 
Lagrangian: 



(3) 



Here, dF^ and are EDM and CEDM, respectively, of a quark flavor /. We perform calculations of the one-loop 
EDM diagrams using mass eigenstates with their effective couplings. The approach frees our numerical results from 
the mass-insertion approximation more commonly adopted in the type of calculations, while analytical discussions 
based of the perturbative diagonalization formulae help to trace the major role of the RPV couplings, especially 
those of the bilinear type. The interesting class of one-loop contributions are obtained from diagrams involving a 
bilinear-trilinear parameter combinations, which were seldom studied by most other authors on RPV physics. The 
bilinear parameters come into play through mass mixings induced among the fermions and scalars (slepton and Higgs 
states), while the trilinear parameter enters a effective coupling vertex. The basic features are the same as those 
reported in the studies of the various related processes [3 Ti. 1391] . Among the latter, our recently available reports on 
6 — !■ s -t- 7 [39} are particularly noteworthy, for comparison. The d quark dipole plays a more important role over that 
of the u quark, when RPV contributions are involved. The b —f s + j diagram is, of course, nothing other than a 
flavor off-diagonal version of a down-sector quark dipole moment diagram. 

Unlike the SM Yukawa couplings, the RPV Yukawa couplings are mostly not so strongly constrained in magnitudes 
[4ll|, and are sources of flavor mixings. A trilinear A'^^^j, coupling couples a quark to another one, of the same or 
different charge, and a generic scalar, neutral or charged accordingly. The coupling together with a SM Yukawa 
coupling at another vertex contributes to the EDMs through some scalar mass eigenstates with RPV mass mixings 
involved, as illustrated in Fig. [T] As the SM Yukawas are flavor diagonal, the charged scalar loops contribute by 
invoking CKM mixings^. Note that under the SVP adopted, the X^jj, parameters have quark flavor indices actually 
defined in the d-sector quark mass eigenstate basis. Analytically, we have the formula for the electric dipole form 
factor for quark flavor /, 



An sin20„ 



M? 



iQf-Qf,)F, 



Ml 

mF 



Qf'F, 



Ml 

mF 



(4) 



For the corresponding flavor off-diagonal transition moment, like b ~> s + ■y , scalar loop contributions do exist [39 
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and for the chromo-electric dipole form factor, 



I 3 



4-iT sin^^ 



m n=l 

Here, Fs{x) and ^4(0;) are Inami-Lim loop functions, given as: 

1 



{-3 + 4x - -2 Inx) 



2(l-a;)3 

^'^""^ " 2(l-a;)3 (l-^' + 2arlna;) . 
For / = -u (i = 1), the coefficients C^;^. are defined by the interaction Lagrangian, 



1-75 , ;?k 1 + 75 



+h.c. , 



(5) 

(6) 
(7) 

(8) 



V* 

1->L* _ ai-n -trin -nl' , jkn -,rik 



92 



92 



yin j)l 



and, for / = d (i = 1), the coefficients C^'J^^ are defined by the interaction Lagrangian, 



1-75 _^ 1+75 



+h.c. , 



(9) 



(10) 



92 



Vdj T/m*7-)i _|_ jki -irnk* 'pfl 
92 92 



X. 



(11) 



with the Y^'^ denoting a sum over (seven) nonzero mass eigenstates of the charged scalar; i.e., the unphysical Goldstone 
mode is dropped from the sum, being the diagonalization matrix, i.e., D^'^M^ = diag{ M^^ , m = 1 -8 } and 
being 8x8 charged-slepton and Higgs mass-squared matrix . Here, i = 1 is the flavor index for the external quark 
while n is for the quark running inside the loop. Note that the unphysical Goldstone mode is dropped from the scalar 
sum because it is rather a part of the gauge; loop contribution, which obviously is real and does not affect the EDM 
calculation. 

In general, there also exist the contributions from neutral scalar loop (involving the mixing of neutral Higgs with 
sneutrino in the loop). These contributions lack the top Yukawa and top mass effects that enhance the contributions 
of charged-slepton loop, and hence less important. The formula for the electric dipole form factor is given as: 



If 



m n=l ^ ' 



(12) 



and for the chromo-electric dipole form- factor is given as. 



/ 3 



m 71=1 



(13) 



where, for / = u (i = 1), the coefficients Kf^'^i are defined by the interaction Lagrangian, 



-C" = 52*(«n)$t(0 



1-75 



1+75 



^(■Ui) +h.c. 



(14) 
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AT'. 



iV: 



6m J J 



(15) 



and, for f = d {i = 1), the coefficients A/"^;^. are defined by the interaction Lagrangian, 



(16) 



52 V2 52 V2 



^(i+2)m - ^'^U+7)rn 



with the denoting a sum over (10) nonzero mass eigenstates of the neutral scalar; i.e., the unphysical Goldstone 
mode is dropped from the sum, D'^ being the diagonalization matrix for the 10 x 10 mass matrix for neutral scalars 
(real and symmetric, written in terms of scalar and pseudo-scalar parts). Again i = 1 is the flavor index for the 
external quark while n is for the quark running inside the loop. 

Having obtained expression for quark EDMs and CEDMs, the next task is to connect them to the hadronic and 
atomic EDMs. This is a non-trivial step owing to non-perturbative effects of QCD for which there is no single 
unambiguous approach. Below we will briefly describe the neutron, deuteron, and mercury EDM formulas that we 
will use for numerical calculations. For details we refer the reader to the excellent review articles 4, 5]. 

For the case of neutron EDM, there have been three different approaches with varying degree of sophistication. The 
simplest of these is the non-relativistic SU{6) quark model [25l. \Wl. [27l. ]. There have also been attempts based 
on the chiral Lagrangian method [sij and QCD sum rule approach [32] ■ In the chiral Lagrangian approach neutron 
EDM is expressed solely in terms of quark CEDM [3l[ : 

d„ = (1.6(l^ + 1.3(f^ + 0.26d;) . (18) 

In the case of QCD sum rule approach one obtains [s^l : 

= (1 ± 0-5) (225MeV)3 ^^■^^^'^d + 0.5<) + 0.7(cf, - 0.25<)] . (19) 

However, the two approaches differ substantially in their dependence on various quark EDMs and and CEDMs and 
their relative importance and hence give very different results. For instance, in the chiral Lagrangian method quark 
EDMs are completely neglected and only CEDMs contribute, while both contribute in the QCD sum rule approach. 
Moreover, strange quark CEDM contribution is dominant in chiral Lagrangian method whereas it is neglected in the 
QCD sum rule approach. Owing to the large differences in these method we use the non-relativistic SU{6) quark 
model and incorporate the quark CEDM contributions using naive dimensional analysis. Although, this approach is 
less sophisticated than the above two methods, it provides a reasonable and reliable enough estimate for our present 
purpose. 

In the SU (6) quark model, one associates a non-relativistic wavefunction to the neutron consists of three constituent 
quarks with two spin states. After evaluating the relevant Clebsch-Gordan coefficients the neutron EDM is written as 

dn ^^{idd-du) , (20) 



where 



d, = rfd-^ + f^dl. (21) 



Here, dF and cF are the contributions to neutron EDM from the electric and chromoelectric dipole operators, respec- 
tively. ?/^(~ 0.61) and rf (~ 3.4) are the respective QCD correction factor from renormalization group evolution 34!|. 
The authors of Ref. i27j| showed that within MSSM chromoelectric dipole form factor contribution to neutron EDM is 
comparable to that from electric dipole form factor and hence should be included. 

Among diamagnetic atoms, mercury atom provides the best limit on atomic EDM, which is a result of T-odd 
nuclear interactions. Schiff screening theorem is violated by finite size effects of the nucleus and are characterized by 
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SchifF moment S which generates T-odd electrostatic potential for atomic electrons. Again the principle approaches 
in the literature are the QCD sum rule and chiral Lagrangian method. The QCD sum rule approach gives |42i | 

dHg^-{dF^-(F^- 0.012 (F^) X 3.2 • 10^^ g _ (22) 

The chiral Lagrangian method gives [i^ 

dHg = -{(fd "(Pu- 0-0051 (F,) X 8.7 X 10-3 e . (23) 

Again, one can see the differences in the two approaches. We shall follow the results of Ref.jil] and interpret the 
bounds on mercury EDM as constraining the combination — cf^| < 2 x lO"^^. As we will soon see in our result 
section, constraints coming from mercury EDM are redundant as neutron and the planned deuteron EDM experiment 
provide much stringent constraints. 

For the case of deuteron EDM, because of rather transparent nuclear dynamics, theoretical uncertainty is much 
smaller. As mentioned in the introduction, the proposed deuteron EDM experiment would reach a sensitivity of 
(1 — 3) X 10"^^ ecm that is about 10 to 100 times better than than the current EDM limits in terms of sensitivity 
to quark EDMs and QCD 9 parameter [l^. Deuteron EDM receives contributions from the constituents proton and 
neutron EDMs as well as CP-odd meson nucleon couplings: 

do = dn + dp + djf " . (24) 

For theoretical predictions again there are two approaches, namely QCD sum rule [is! ] and the chiral Lagrangian 
approach [i^. Authors of Ref.ji^ show that in the SU{2) chiral perturbation theory proton and neutron EDMs 
exactly cancel each other and hence compute the EDM using QCD sum rule approach. Authors of Ref.[4G] have 
shown that this is not true after introducing strange quark and corresponding meson. However, the two approaches 
agree to within error-bars regarding the dominant contributions due to quark CEDMs and following Ref. [33| we take 
deuteron EDM to be doidq, d^) ~ —e {d^ — d^) 5^3^ and use the best-fit value for our analysis. 



IV. RESULTS AND DISCUSSIONS 



Prom the discussion in previous section it is clear that analysis of neutron, deuteron, and mercury EDMs boils 
down to understanding quark EDMs and CEDMs. Quark EDMs involve violation of CP but not R-parity. Thus RPV 
parameters should come in combinations that conserve R-parity. From an inspection of formulae it is clear that the 
contributions from two A' couplings cannot lead to EDM as these violate lepton number by two units which has to 
be compensated by Majorana like mass insertions for neutrino or sneutrino propagators ^. If one is willing to admit 
more than two A' to be non-zero, than in principle one can have contributions to EDM at one loop level from the 
Majorana like mass insertions but these would be highly suppressed. With only two RPV couplings, the only other 
possibility at one loop level is to have a combination of trilinear and a bilinear (/i^, Bi or m'^^.) RPV couplings in 
such a way that lepton number is conserved. However not all the three bilinears mentioned above are independent as 
they are related by tadpole relation in the single VEV paramctrization [2l[ 

Bi tan/3 — mlg. + fJ-lfJ-i ■ (25) 

Using the above tadpole equation we eliminate "t-^,,. in favor of /ii and Bi. Contributions from the combination 
l^iKjk^ through squark loops, had been extensively studied in ^| in detail. Here, we shall focus on the B*\[-f, kind 
of combination as illustrated in FiglT] Such a combination can contribute through charged-scalar (charged-slepton, 
charged Higgs mixing) loop or neutral scalar (sneutrino, neutral Higgs mixing) loop. The fermions running inside the 
loops are quarks, instead of the gluon or a colorless fermion as in the case of the squark loops. Before we discuss the 
numerical results it would be worthwhile to discuss the analytical results which can then be compared with numerical 
results. _ _ 

Let us focus on RPV part of Im(C^^^.C^^;), in Eq.(l4|) for the case of d quark EDM. It is given as 

Im(Ct„,C:J„P.=Im[(2/„„y,'^i„I?l:j X {\',kiV-^*V\^^2)j] ■ (26) 



^ A combination A'*, A' ., is real and hence does not contribute to EDM 
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For the u quark dipole, we have 



(27) 



I* 



Interestingly, the entries of the slepton-Higgs diagonahzing matrix involved are the same in both terms above. Bilinear 
RPV terms are hidden inside the slepton-Higgs diagonalization matrix elements. To see the explicit dependence on bi- 
linear RPV terms let us look at the diagonahzing matrix elements of charged-slepton Higgs mass matrix, T^'(^j^2)m'^ 
more closely. When summed over the index m, it of course gives zero, owing to unitarity. This is possible only for the 
case of exact mass degeneracy among the scalars when loop functions factor out in the sum over m scalars. In fact, 
the final result involves a double summation over the m scalar and n fermion (quark) mass eigenstates. Either the 
unitarity constraint over the former sum or the GIM cancellation over the latter predicts a null result whenever the 



mass dependent loop functions F4 



M. 



I'n. 



M? 



and F-? 



M, 



f'r. 



M? 



can be factored out of the corresponding summation due to 



mass degeneracy. In reality however, these elements are multiplied by non-universal loop functions giving non-zero 



EDM. Restricting to first order in perturbation expansion of the mass-eigenstates, 1^\j+2)m'^im non-zero for m = 2 



and m — j + 2, both giving similar dependence on RPV but with opposite sign 



state that is dropped from the sum here). Also V 



I 

(i+2),2 



^l(i+2)- 



1 is the unphysical Goldstone 



For m — 2 one obtains [211 



V 



I 

0+2)2 



^xO(l) 



(28) 



Here, AI^ denotes the difference in the relevant diagonal entries (generic mass-squared parameter of the order of soft 
mass scale) in charged-slepton Higgs mass matrix. To first order in perturbation expansion there is no contribution to 
EDM from fii . If one goes to second order in perturbation expansion, one gets a contribution from the term m = j + 5 
which is given as [2l[ 



V' 



I* 



0+2)0+5) ^10+5) 



{Al - Ho tan P)mj V2M„ sin /3{nkKjj ) 



52M2 



(29) 



There are a few important things to be noted here 
1 



Notice that /x^ enter at second order in perturbation expansion. Moreover they are accompanied by corresponding 
charged- lepton mass-squared and hence /i*A^^ j, contributions are always suppressed (later below we will elaborate 
on this). Also notice that in principle trilinear parameter Xkjj also contribute through fJ-k^kjj- But they have 
to be present in addition to the trilinear parameter X'^Jj^, thus making it a fourth-order effect in perturbation 
and hence negligible. Thus, we will focus on B*X'-j^ effects which are interesting. 

2. Even if all RPV parameters are real, CKM phase in conjunction with real RPV parameters could still induce 
EDM. As we will soon see, this could be sizable. 

3. It is clear that d quark EDM receives much larger contribution owing to top- Yukawa and proportionality to top 
mass. There are nine trilinear RPV couplings X'--^ that contribute to d quark EDM. X'^^^ has the largest impact 



owing to the least CKM suppression. 



The nit enhancement feature is not there in the two loop contributions 



4. 



5. 



All the 27 trilinear RPV couplings (A^^,.) contribute to u quark EDM. However, the absence of enhancement 
from the top mass in the loop and the uniform proportionality to up- Yukawa considerably weakens the type of 
contribution to u quark EDM. 

There is no RPV neutral scalar loop contribution to u quark EDM. However, there are one-loop contributions 
to d quark EDM from the neutral-Higgs sneutrino mixing due to the combination of BiX'^n with Majorana like 
mass-insertion in the loop^. From the EDM formula, it is clear that this would be about the same magnitude 
as the u quark EDM due to charged-Higgs charged slepton mixing and hence highly suppressed. 



^ The Majorana like mass- insertion can be considered a result of the non-zero Bi. It manifest itself in our exact mass eigenstate 
calculation as a mismatch between the corresponding scalar and pseudo-scalar parts complex "sneutrino" state which would otherwise 
have contributions canceling among themselves. With the non-zero Bi, the involved diagram is one with two X'^n coupling vertices and 
a internal d quark. Hence, the type of contribution is possible only with the single A' coupling. 
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From the above analytical discussion, we illustrated clearly how various combinations of trilinear and bilinear RPV 
parameters contribute to neutron EDM. We now focus on the numerical results. In order to focus on individual 
contributions we keep a pair of RPV coupling (a bilinear and a trilinear) to be non-zero at a time to study its impact. 
We have chosen all sleptons and down-type Higgs to be 100 GeV (up-type Higgs mass and Bq being determined from 
electroweak symmetry breaking conditions), /Xq parameter to be -300 GeV, A parameter at the value of 300 GeV, 
and tan/3 = 3 (we will show the impact of some parameter variations below). Since the A' couplings are on the same 
footing as the standard Yukawa couplings, they are in general complex. We admit a phase of tt/A for non-zero A' 
couplings while still keeping the CKM phase. The phases of the other R-parity conserving parameters are switched 
off. Since it is the relative phase of the B*\^-f, product that is important, we put the phase for Bi to be zero without 
loss of generality. We do not assume any hierarchy in the slcptonic spectrum. Hence, it is immaterial which of the 
bilinear parameter Bi is chosen to be non-zero. We choose B^ to be non-zero but all the results hold good for Bi and 
B2 as well. 

In Fig. ([21) we have plotted the neutron, deuteron, and mercury EDMs, as well as d-quark EDM and CEDM, against 
the most important combination Im^B'^X'^^i) normalized by /i^. This combination has the largest impact owing to top 
Yukawa and top mass dependence. We have not plotted u-quark EDM as it is about 5 orders of magnitude smaller 
than d-quark EDM (reflecting up-top hierarchy). The description of various lines and symbols is given in the caption 
for the figure. It is straight forward to understand the features of the plot. From the Eqs.Q and ^ given that 
loop functions F3 w F4 for the parameters considered here, one can see that d^^/dr^ ~ 5/(2 g^) « 2.5. This can be 
roughly observed in the plot. The plot also shows that total neutron EDM is same as the d-quark EDM contribution. 
To understand this, observe that u-quark EDM is negligible and hence the total neutron EDM is sum of d-quark 
EDM and CEDM. CEDM suffers not only a factor of 2.5 suppression as discussed above, but also has additional e/An 
suppression, which more than compensates QCD enhancement factor of 3.4. For d-quark EDM contribution, QCD 
suppression factor is compensated by the Clebsch-Gordan coefficient factor of 4/3 and the total neutron EDM works 
out to be roughly same as d-quark EDM. In the similar way one can see from the formulas for deuteron and mercury 
EDMs that dng ~ c^/3.14 and ddcutcron = 1-6 x d^. All these relations are born out in the plot. 

Coming to bounds, from Figl^l we see that the null result at the proposed Los Alamos EDM experiment can lead 
to a particularly stringent constraint on Im(_B3 A332)//io < 1.4 x 10^^. For the same combination the bound we obtain 
from the sensitivity of the proposed deuteron EDM experiment is 4.0 x 10~^ and from the present best limits on 
neutron EDM and Hg EDM experiment, 4.0 x 10^ and 1.9 x 10~^ respectively. This gives us an idea about the 
unprecedented sensitivity of the proposed Los Alamos neutron EDM experiment and the deuteron EDM experiment. 

FigO shows the contours of neutron EDM in the plane of magnitudes for the couplings X':jji and B3 with the relative 
phase fixed at 7r/4. The contour for present experimental bound is shown in pink dotted line. The plot shows that a 
sizable region of the parameter space are ruled out. Successive contours show smaller values of neutron EDM with the 
smallest being 10~^® e cm, the projected sensitivity of Los Alamos experiment. FigU] shows the contours of deuteron 
EDM in the plane of magnitudes for the couplings X'^^i and B-s with the relative phase fixed at 7r/4. Region between 
pink dotted line and green dashed line is the expected sensitivity of the proposed deuteron experiment. Black dashed 
line shows a possible order of magnitude improvement in this. One can see that any null result can rule out huge 
regions in parameter space. 

So far we have kept certain parameters like the phase of the RPV combination, the fj,Q parameter and the slepton 
mass spectrum fixed. To get a better understanding of the allowed regions in the overall parameter space, let us 
focus on variations of the parameters one at a time. In Fig. ([5]) we have plotted neutron EDM versus the slepton 
mass parameter niL = ttie (with I-B3I = 200 GeV^, IA332I = .05 and relative phase of 7r/4). As expected the neutron 
EDM falls with the increasing slepton mass. More careful checking reveals that the result is sensitive only to one 
slepton mass parameter, the mass of the third L-handed slepton here. It is also easy to understand from our analytical 
formulas that the dominating contributions among the various scalar mass eigenstates for the case of B*\'^-j, come 
from the j-th L-handcd slepton and the Higgs. Various horizontal lines arc the experimental bounds. It is interesting 
to note that even if slepton mass is as high as 1 TeV, it still gives neutron EDM contribution that is an order of 
magnitude larger than the Los Alamos experiment sensitivity. In Fig. ([6]) we have shown the variation of neutron 
EDM with the /Xq parameter. Although parameter /Xq does not directly figure in the EDM formula, its infiuence is 
felt through the Higgs spectrum. Larger /io leads to heavier Higgs spectrum which suppresses the EDM contribution. 
Again, it is interesting to see that /xo = ITeV, still gives neutron EDM contribution that is larger than the sensitivity 
of Los Alamos experiment. The large value of neutron EDM (much larger than the reach of Los Alamos experiment) 
even for slepton mass and /io parameter at 1 TeV is not only encouraging but also indicates that if a non-zero EDM 
is measured, one should exercise caution in the interpretation in terms of any new physics model. Our analytical 
formulas show that there is no strong dependence on tan/3 (we have checked this numerically as well) and hence 
we have kept tan/3 — 3 fixed in all the plots. Relative insensitivity to the value of tan/3 here can be contrasted to 
certain interesting predictions like Branching fraction for Bs — > /x^/i^ in MSSM which is boosted by three orders 
of magnitude in large tan/3 region. Any constraints on tan/3 following non-observations of Bs — > fJ-^fJ-^ would still 
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TABLE I: Here we list the normalized upper bounds for several combinations of bilinear and trilinear RPV parameters, with 
some variation in the input parameters. The bounds essentially depend on the values of parameters like rrii, /xq and rrijf^ (tih^ 
and Bo being determined from EW symmetry breaking condition), tan/3 has been kept fixed at 3 as EDM has a very mild 
dependence on tan/3. Note that all these bounds include CEDM contributions. Columns headings I,II,III refer to bounds from 
current neutron EDM experiment, from the possible null results at the Los Alamos neutron EDM experiment, and the possible 
null results at the deuteron EDM experiment respectively. As is clearly seen, any null results at Los Alamos experiment can 
lead to bounds that are over two orders of magnitude larger than the current bounds 
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"For this case the constraints correspond to the real part of RPV combination as the relative phase is zero. 



leave our predictions unaffected. The fact that our predictions largely depend on hitherto unconstrained combination 
B*\'^ji and already known top Yukawa and top mass, makes it all the more interesting. 

To go beyond the illustrative case of the parameter combination B*\^^^ , we list in table [T] the bounds due to the 
neutron EDM from the current Grenoble results (in the column under I), due to plausible null result at Los Alamos 
experiment (in the column under II), as well as due to any null results of the future deuteron EDM experiment (in 
the column under III), on the nine combinations |i3,*A^j]^| normalized by (100 GeV)^. It is interesting to note the 
difference in bounds for case (a) and (b) for coupling combination |i?*A^3]^| and |i?*A^2il- The inputs for Case (b) 
are identical to case (a) except for RPV phase being zero for case (b). Case (b) thus relies solely on CKM phase. 
Interestingly the bound for |i3*A-3;^| changes very marginally from case (a) to (b) whereas the bound for |i?*A-2il 
weakens by about an order of magnitude. To understand this difference in behavior of |Sj*A^3]^| and |i3i*A^2ili with 
and without a complex phase in the RPV couplings, we have plotted in Fig. [7] the allowed region in the plane of 
relative phase of S3A331 and the |i?3A33]^| (left) and relative phase of B3A321 and the IB3A321I (right). One can see 
that the bound for |S3A33]^| is about the same for relative phase in B'i^X'^^^ of or 7r/4 but the bound for IS3A321I 
strengthens by about an order of magnitude as the phase in the RPV coupling increases from to 7r/4 suggesting a 
collaborative effect between the CKM phase and the RPV phase. Fig[8] is similar to FiglT] except that it is now for 
deuteron EDM. We see that a similar pattern here also because of the similar qualitative dependence on down-quark 
EDM, implying similar interference pattern between CKM phase and RPV phase induced contributions. For the case 
of i3*A^2i there is no CKM phase involved. In the table we have fixed the sign of /i,, parameter to be negative. In 
the Figl6] it is seen that absolute value of neutron EDM is symmetric with respect to sign of /Iq parameter and hence 
positive /io should give identical bounds. The variation of bounds in table |T] with changes in parameters /.to and the 
slepton and Higgs mass very much follows the pattern found in plots discussed earlier. 

Before we conclude, we would like to briefly comment on two things. In the table HI we have only mentioned 
the bounds for the nine combinations |i?*A^j]^|, whereas earlier in the text we did mention that in principle all the 
twenty-seven X'^-^, together with bilinear couplings do contribute to the neutron EDM. The couplings other than X^^-^ 
contribute to u quark EDM and all the contributions are proportional to up- Yukawa (in contrast to the presence 
of a contribution proportional to top- Yukawa for the d quark EDM). Strength of the corresponding contributions is 
substantially weaker (typically by 5 to 6 orders of magnitude) and hence do not lead to meaningful bounds. The 
other thing is about the possible fi* ■ X[j^ contributions. It can be seen in Eq. ([29|) that these are second order in 
perturbation and also accompanied by lepton mass jii. Thus, a typical /x*A^^^ contribution is substantially weaker 
than the corresponding B*X^jf, contribution. For instance, with the similar inputs for other SUSY parameters as 
described above, if one takes fi^ — 10^^ GeV (dictated by requirement of sub-ev neutrino masses), A33;^ — .05 and a 
relative phase of 7r/4, one obtains neutron EDM of 6.0 • 10"'^^ which is about six orders of magnitude smaller than the 
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present experimental bound on neutron EDM. If one goes by the upper bound on the mass of Vr of 18.2 MeV \^f\ , 
fis could be as large as 7 GeV for tan/3 = 2 and sparticle mass ^ 300 GeV 20;]. For fj,3 = 1 GeV we obtain neutron 
EDM of 6.1 • 10-2^, still about three orders of magnitude smaller than the experimental bound. These numbers can 
be easily compared with the fJ-i^'m contributions to the d quark EDM through chargino loop in Table 1 of ref . [37j . 
There the corresponding contribution (with /ig = 1 GeV) is much weaker (of order 10~^^) as it lacks the top Yukawa 
and top mass enhancements. In the same table one finds that corresponding contribution to gluino loop is much 
stronger (of order 10~^^) due to proportionality to gluino mass and strong coupling constant. In the light of above 
reasons one can appreciate that the contributions due to soft parameter Bi are far more dominating in the present 
scenario of quark loop contributions. 

V. CONCLUSIONS 

We have made a systematic study of the influence of the combination of bilinear and trilinear RPV couplings on 
the neutron, deuteron, and mercury EDMs, including the contributions from chromo-electric dipole form factor. Such 
combinations are interesting because this is the only way RPV parameters contribute at one-loop level to EDM. The 
fact that the form factors are proportional to top Yukawa and top mass with the hitherto unconstrained combination 
of RPV parameters, makes them all the more interesting. Such class of diagrams have a quark as the fermion running 
inside the loop. EDMs of neutron, deuteron, and mercury are ultimately expressed in terms of quark EDMs and 
CEDMs. In our analytical expressions obtained based on perturbative diagonalization of the scalar mass-squared 
matrices, we demonstrated that charged-slepton Higgs mixing loop contribution to d quark EDM far dominates the 
other contributions due to a diagram with the top-quark in the loop. In our numerical exercise we have obtained 

robust bounds on the combinations ^(^Q^^Q^Yy ~ 1:2,3 that have not been reported before. The bounds 

are reported for the current best limits on neutron EDM as well as projected for any null results in the future 
improvements at Grenoble and Los Alamos experiments as well as deuteron EDM experiment. It turns out that 
measurements of mercury EDM are not as strongly constraining as current neutron EDM limits, whereas any null 
result at deuteron EDM and/or Los Alamos neutron EDM experiment can lead to much stronger constraints on RPV 
parameter space (for instanccIm(S3A33]^)//iQ < 1.4 x 10~^). Even if the RPV couplings are real, they could still 
contribute to quark EDMs via CKM phase. For some cases CKM phase induced contribution is as strong as that due 
to an explicit complex phase in the RPV couplings. We find that even if slepton mass or fig are as high as 1 TeV, it 
could still lead to bounds that are well within the reach of Los Alamos neutron EDM experiment. There also exist 
contributions involving /i*A^jj,. However these are higher order effects which are further suppressed by proportionality 
to charged lepton mass. Since /i^ are expected to be very small (of order 10"'^ GeV) for sub-eV neutrino masses, such 
contributions are highly suppressed. Our results presented here make available a new set of interesting bounds on 
combinations of RPV parameters. 

One further thing worth some attention here is the question of to what extent our choice of model formulas for the 
hadronic EDMs from quark dipoles, as discussed and justified above, affect our major conclusions One can think 
about presenting full numerical results and comparison for the cases of the various different model formulas. However, 
we find that neither feasible nor desirable. We did do some numerical checking, but consider a brief summary here 
as the only appropriate thing to present. Two special features of the RPV model dictate the answer. The hadron 
EDMs are to be determined from the EDMs and CEDMs of the u, d, and s quark. For the RPV model, the resulted 
dipoles for w-quark are much smaller than of the d-quark. For the s-quark, however, dominant contributions involve 
similar diagrams with trilinear RPV couplings of different familiy indices, typically replacing a 1 for the d by a 2 for 
the s. The latter complicaton makes it difficult to fully address the impact within the scope of the present study. 
However, one can at least see that for the RPV parameter combinations playing a dominant role in generating d 
dipoles, our main focus here, their contributions to the s dipoles are going to be unimportant. Neglecting the the 
s and u dipoles, models formulas from QCD sum rule and the valence quark model give roughly c?„ ~ 0.92 d^ and 
dn — 0.96 d^, respectively. The case for comparison with the chiral lagrangian approach is more complicated. But the 
result agrees within a factor of three with the two numbers. The deuteron EDM numbers from the two approaches 
agree within error bars. And the Mercury EDM constraint is always negligible when compared to that of neutron 
and Deuteron, even when s dipoles are to be considered. Hence, our results here are not much affected by the hadron 



* The question is prompted by a journal referee, to whom we express our gratitute. 
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FIG. 1: d quark EDM due to _B*Ai3i combination, 
contribution. 



Due to top Yukawa and top-mass dependence this is the most dominant 
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FIG. 2: (color online) The neutron EDM (in units of e.cm) versus the combination ^(Bl A331) normalized by /ig.The horizontal 
lines show present EDM bounds as well as improvements possible in the future experiments. red solid line is current experimental 
bound for neutron EDM 13], pink dotted line is current experimental bound for mercury EDM [l3|, green dashed line is the 
goal expected to be reached by the current Grenoble neutron EDM experiment [T^ . cyan dash-dotted line is the expected 
measurement of of the proposed deuteron EDM experiment [l^, and blue hashed line is the expected improvement at the 
proposed Los Alamos neutron EDM experiment [l^ . The red filled squares stand for the total neutron EDM prediction, green 
filled circles are d-quark EDM, blue open circles are d-quark CEDM, cyan filled squares are predictions for deuteron EDM and 
magenta filled triangles are predictions for mercury EDM, in the present framework described in the previous section. 



14 




\X 



3311 



5e-05 0.0001 0.00015 0.0002 0.00025 



IBgl/n'o 



FIG. 3: A contour plot for various values of neutron EDM in e.cm in the plane of real A331 and B^. (with a relative phase of 
7r/4). Cyan dash-dotted line is the previous best experimental bound, pink dotted line is the current best bound, blue hashed 
line is the expected reach of Grenoble experiment, black dashed line is the Los Alamos projected sensitivity. One can see a 
sizable improvement in the regions ruled out in parameter space with improvement in sensitivity. 




FIG. 4: A contour plot for various values of deuteron EDM in e.cm in the plane of real A331 and B3 (with a relative phase of 
7r/4). Region between pink dotted line and green dashed line is the expected sensitivity of the proposed deuteron experiment. 
Black dashed line shows a possible order of magnitude improvement in this. Again one can see that any null result can rule 
out huge regions in parameter space. 
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FIG. 5: Various EDMs vs. the slepton mass parameter ttil ~ ms (with IBs] — 200 GeV^, IA331I — .05 and relative phase of 
7r/4). Horizontal lines are the experimental bounds. 
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FIG. 6: Various EDMs vs. fj,o parameter with \Bs \ = 200 GeV'^, |A33i[ = .05 and relative phase 7r/4. Horizontal lines are the 
experimental bounds. 
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FIG. 7; Dotted region is the parameter space allowed by neutron EDM upper bound in the plane of relative phase of 
-^331 -B3 and the |A33ii?|| (top) and relative phase of A321S3 and the IA321-B3I (bottom). 



17 



CM 



CO 

CD 

CO 
CO 



1 .2e-05 
1e-05 
8e-06 
6e-06 
4e-06 
2e-06 




T — I 1 — I 1 — I 1 — I 1 — I — I 1 — I 1 — I 1 — I 1 — I — I 1 — I 1 — I 1 — I 1 — I — I 1 — I 1 — I 1 — T" 




■150 -100 -50 50 100 150 

5 

Relative phase of ^331 B3 



CO 

00 



CO 
CO 



■''u 

v:. 



'■>>• v. 



6e-05 
5e-05 
4e-05 
3e-05 
^ 2e-05 
1e-05 


-150 -100 -50 50 100 150 
Relative phase of ^331 B3 

FIG. 8: Dotted region is the parameter space allowed by proposed deuteron EDM upper bound in the plane of relative 
phase of A331BI and the IA331BII (top) and relative phase of A321-B3 and the IA321B3I (bottom). 
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